On the final spin from the coalescence of two black holes 
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We provide a compact analytic formula to compute the spin of the black hole produced by the coalescence 
of two black holes following a quasi-circular inspiral. Without additional fits than those already available for 
binaries with aligned or antialigned spins, but with a minimal set of assumptions, we derive an expression that 
can model generic initial spin configurations and mass ratios, thus covering all of the 7-dimensional space of 
parameters. A comparison with simulations already shows very accurate agreements with all of the numerical 
data available to date, but we also suggest a number of ways in which our predictions can be further improved. 
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The evolution of black hole binary systems is one of the 
most important problems for general relativity, and more re- 
cently for astrophysics, as such systems enter the realm of ob- 
servation. Recent advances in numerical relativity have made 
it possible to cover the entire range of the inspiral process, 
from large separations at which post-Newtonian (PN) calcula- 
tions provide accurate orbital parameters, through the highly 
relativistic merger, to ringdown. For many studies of astro- 
physical interest, such as many-body studies of galactic merg- 
ers, or heirarchical models of black-hole formation however, 
it is impractical to carry out evolutions with the full Einstein, 
or even post-Newtonian, equations. Fortunately, recent binary 
black-hole evolutions in full general relativity have shown that 
certain physical quantities can be estimated to good accuracy 
if the initial encounter parameters are known. In particular, 
this paper develops a rather simple and robust formula for de- 
termining the spin of the black-hole remnant resulting from 
the merger of rather generic initial binary configurations. 

To appreciate the spirit of our approach it can be conve- 
nient to think of the inspiral and merger of two black holes 
as a mechanism which takes, as input, two black holes of ini- 
tial masses Mi, Mi and spin vectors Si, S2 and produces, as 
output, a third black hole of mass Mg n and spin S& n . In con- 
ditions of particular astrophysical interest, the inspiral takes 
place through quasi-circular orbits since the eccentricity is 
removed quickly by the gravitational -radiation reaction 
Furthermore, at least for nonspinning equal-mass black holes, 
the final spin does not depend on the value of the eccentricity 
as long as it is not too large The determination of Mg n 
and Sfl n from the knowledge of M12 and Si^, is of great 
importance in several fields. In astrophysics, it provides infor- 
mation on the properties of isolated stellar-mass black holes 
produced at the end of the evolution of a binary system of 
massive stars. In cosmology, it can be used to model the dis- 
tribution of masses and spins of the supermassive black holes 
produced through the merger of galaxies (see ref. ^ for an 
interesting example). In addition, in gravitational-wave as- 
tronomy, the a-priori knowledge of the final spin can help the 
detection of the ringdown. What makes this a difficult prob- 
lem is clear: for binaries in quasi-circular orbits the space of 
initial parameters for the final spin has seven dimensions (i.e., 



the mass-ratio q = M2/M1 and the six components of the spin 
vectors). A number of analytical approaches have been devel- 
oped over the years to determine the final spin, either exploit- 
ing the dynamics of point-particles 10, or the PN approxi- 
mation [6], or using more sophisticated approaches such as the 
effective-one -body approximation J3]. Ultimately, however, 
computing ag n = Sa n /Ml n accurately requires the solution 
of the full Einstein equations and thus the use of numerical- 
relativity simulations. Several groups have investigated this 
problem over the last couple of years MI^ [Tfl[Tll[l2ilT3ll . 

While the recent possibility of measuring accurately the fi- 
nal spin through numerical-relativity calculations represents 
an enormous progress, the complete coverage of the full pa- 
rameter space uniquely through simulations is not a viable 
option. As a consequence, work has been done to derive 
analytic expressions for the final spin which would model 
the numerical-relativity data but also exploit as much infor- 
mation as possible either from perturbative studies, or from 
the symmetries of the system 11 [H [H [H [H. In 
this sense, these approaches do not amount to a blind fitting 
of the numerical-relativity data, but, rather, use the data to 
construct a physically consistent and mathematically accurate 
modelling of the final spin. Despite a concentrated effort in 
this direction, the analytic expressions for the final spin could, 
at most, cover 3 of the 7 dimensions of the space of parame- 
ters lfl3ll . Here, we show that without additional fits and with 
a minimal set of assumptions it is possible to obtain the ex- 
tension to the complete space of parameters and reproduce all 
of the available numerical-relativity data. Although our treat- 
ment is intrinsically approximate, we also suggest how it can 
be improved. 

Analytic fitting expressions for a^ n have so far been built 
using binaries having spins that are either aligned or an- 
tialigned with the initial orbital angular momentum. This is 
because in this case both the initial and final spins can be pro- 
jected in the direction of the orbital angular momentum and 
it is possible to deal simply with the (pseudo)-scalar quanti- 
ties a\, a,2 and a^ n ranging between —1 and +1. If the black 
holes have equal mass but unequal spins that are either par- 
allel or antiparallel, then the spin of the final black hole has 
been shown to be accurately described by the simple analytic 
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fit 

afin(ai,a 2 ) = Pa +pi(ai + a 2 ) +p 2 (ai + a 2 ) 2 , (1) 

where p = 0.6883 ± 0.0003, pi = 0.1530 ± 0.0004, and 
p 2 = —0.0088 ± 0.0005. When seen as a power series of the 
initial spins, expression (Q~|) suggests an interesting physical 
interpretation. Its zeroth-order term, in fact, can be associated 
with the (dimensionless) orbital angular momentum not radi- 
ated in gravitational waves and amounting to ~ 70% of the 
final spin at most. The first-order term, on the other hand, can 
be seen as the contributions from the initial spins and from 
the spin-orbit coupling, amounting to ~ 30% at most. Finally, 
the second-order term, includes the spin-spin coupling, with a 
contribution to the final spin which is of ~ 4% at most. 

If the black holes have unequal mass but spins that are equal 
and parallel, the final spin is instead given by the analytic 

fit 01 



dfin(o, v) = a + s^o^v + s^av 2 + t^av + 
2V3v + t 2 v 2 + t 3 i> 3 , 



(2) 



where v is the symmetric mass ratio v = M\Mij (Mi+M 2 ) 2 , 
and where the coefficients take the values S4 = — 0.129 ± 
0.012, s 5 = -0.384 ± 0.261, to = -2.686 ± 0.065, t 2 = 
-3.454 ± 0.132, U = 2.353 ± 0.548. Although obtained 
independently in flUl and lfl3ll . expressions (fl3 and (O are 
compatible as can be seen by considering (f2]i for equal-mass 
binaries (v = 1/4) and verifying that the following relations 
hold within the computed error-bars 
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As long as the initial spins are aligned (or antialigned) with 
the orbital angular momentum, expression (O can be extended 
to unequal-spin, unequal-mass binaries through the substitu- 
tion 
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To obtain this result, it is sufficient to consider ([TJ and (f2| as 
polynomial expressions of the generic quantity 
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where a to t = (ai + fl2<7 2 )/ (1 + <?) 2 is the total dimensionless 
spin for generic aligned binaries. In this way, expressions ([TJ 
and (f2]l are naturally compatible, since a = (ai + a 2 )/2 for 
equal-mass unequal-spin binaries, and a = a for unequal- 
mass equal-spin binaries. Furthermore, the extreme mass- 
ratio limit (EMRL) of expression (f2]i with the substitution 
yields the expected result: afi n (ai, a 2 , v = 0) = a\. 

As already commented above, the predictions of expres- 
sions d2J and cover 3 of the 7 dimensions of the space of 
parameters for binaries in quasi-circular orbits; we next show 
how to to cover the remaining 4 dimensions and derive an an- 
alytic expression for the dimensionless spin vector a^ n of the 



black hole produced by the coalescence of two generic black 
holes in terms of the mass ratio q and of the initial dimension- 
less spin vectors ai, 2 . To make the problem tractable analyti- 
cally, 4 assumptions are needed. While some of these are very 
natural, others can be relaxed if additional accuracy in the es- 
timate of afi n is necessary. It should be noted, however, that 
removing any of these assumptions inevitably complicates the 
picture, introducing additional dimensions, such as the initial 
separation in the binary or the radiated mass, in the space of 
parameters. 

As a result, in the simplest and yet accurate description the 
required assumptions are as follows: 

( i) The mass radiated to gravitational waves Af ra( j can be 
neglected i.e., M fin = M = M-y + M 2 . We note that 
M rad /M = 1 - Mfi n /M « 5 - 7 x 10~ 2 for most of the 
binaries evolved numeri cally . The same assumption was ap- 
plied in the analyses of ll ll [l3ll . as well as in Q], Relaxing 
this assumption would introduce a dependence on Mfi n which 
can only be measured through a numerical simulation. 

(ii) At a sufficiently large but finite initial separation the 
final spin vector can be well approximated as the sum of 
the two initial spin vectors and of a third vector i 



S'fin — S\ + S 2 



(6) 



Differently from refs. |01 and JH], where a definition similar to 
(|6]l was also introduced, here we will constrain £ by exploiting 
the results of numerical-relativity calculations rather than by 
relating it to the orbital angular momentum of a test particle 
at the innermost stable circular orbit (ISCO). When viewed as 
expressing the conservation of the total angular momentum, 
eq. (O also defines the vector £ as the difference between the 
orbital angular momentum when the binary is widely sepa- 
rated L, and the angular momentum radiated until the merger 

J rad, IS., £ = L J ra d. 

(Hi) The vector £ is parallel to L. This assumption is cor- 
rect when Si = — S 2 and q = 1 [this can be seen from the PN 
equations at 2.5 order], or by equatorial symmetry when the 
spins are aligned with L or when Si = S 2 = (also these 
cases can be seen from the PN equations). For more general 
configurations one expects that £ will also have a component 
orthogonal to L as a result, for instance, of spin-orbit or spin- 
spin couplings, which will produce in general a precession of 
£. In practice, the component of £ orthogonal to L will cor- 
respond to the angular momentum J^ ad radiated in a plane 
orthogonal to L, with a resulting error in the estimate of \£\ 

which is - |J r X ad | 2 /W 2 - \Ji d \ 2 /(2V3M 1 M 2 ) 2 M. Al- 
though these errors are small in all the configurations that we 
have analysed, they may be larger in general configurations. 
Measuring J rad via numerical-relativity simulations, or esti- 
mating it via high-order PN equations, is an obvious way to 
improve our approach. A similar assumption was also made 
inref. @|. 

(iv) When the initial spin vectors are equal and opposite 
(Si = —S 2 ) and the masses are equal (q = I), the spin of 
the final black hole is the same as for the nonspinning bina- 
ries. Stated differently, equal-mass binaries with equal and 
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opposite-spins behave as nonspinning binaries, at least when 
it comes down to the properties of the final black hole. While 
this result cannot be derived from first principles, it reflects 
the expectation that if the spins are the same and opposite, 
their contributions to the final spin cancel for equal-mass bi- 
naries. Besides being physically reasonable, this expectation 
is met by all of the simulations performed to date, both for 
spins aligned with L ifTH [l3ll and orthogonal to L ifioll . In 
addition, this expectation is met by the leading-order contri- 
butions to the spin-orbit and spin-spin point-particle Hamil- 
tonians and spin-induced radiation flux 0, [Till . A similar 
assumption is also made, although not explicitly, in ref. JH] 
which, for Stot = 0, predicts l = and |afi n | = £ rb(t = 
0, |ofi n |)/M = const, [cf. eqs. (12)-(13) in ref. [0]. 

Using these assumptions we can now derive the analytic 
expression for the final spin. We start by expressing the vector 
relation (O as 
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(ai + a 2 q 2 + £q) , 



(7) 



where a fin = S^ n /M 2 [cf. assumption (i)], £ = £/(MiM 2 ), 
ai,2 = Si, 2 /M 2 2 , and its norm is then given by 



afi n 



(1 + 9) 2 



a 2 |V + 2|a 2 ||ai|<f cosa- 



2 (|oi| cos/3 + \a 2 \q 2 cos 7 ) \£\q + \£\ 2 q 2 



1/2 



(8) 



where the three (cosine) angles a, j3 and 7 are defined by 

cos a = &i ■ a 2 , cos f3 = ■ £ , cos 7 = a 2 ■ £ . (9) 

Because ai 2 || Si t 2 and t || L [cf. assumption (Hi)], the an- 
gles a, /3 and 7 are also those between the initial spin vectors 
and the initial orbital angular momentum, so that it is pos- 
sible to replace 01,2 with §i. 2 and £ with L in (O. Note 
that a, /3 and 7 are well-defined if the initial separation of the 
two black holes is sufficiently large [cf. assumption ( ii)] and 
that the error introduced by assumption ( Hi) in the measure of 
cos a, cos /3 and cos 7 is also of the order of | J^ d |/|£|. 

The angle £>fi n between the final spin vector and the ini- 
tial orbital angular momentum can be easily calculated from 
I afi n I . Because of assumption (Hi), the component of the final 
spin in the direction of L is [cf. eq. ©] 



a fin = "An ' ^ 



ai|cos/?+ |a 2 |g 2 cos7+ \£\q ^ 



(1 + 9) 



so that cos 6> nn = a nn / 1 \ , and the component orthogonal to 
the initial orbital angular momentum is a^ n = |afi n | sin #fi n . 

In essence, therefore, our approach consists of considering 
the dimensionless spin vector of the final black hole as the 
sum of the two initial spins and of a third vector parallel to the 
initial orbital angular momentum when the binaries are widely 
separated. Implicit in the assumptions made, and in the logic 
of mapping an initial-state of the binary into a final one, is the 
expectation that the length of this vector is an intrinsic "prop- 
erty" of the binary, depending on the initial spin vectors and 



mass ratio, but not on the initial separation. This is indeed 
a consequence of assumption (ii): because the vector £ mea- 
sures the orbital angular momentum that cannot be radiated, 
it can be thought of as the angular momentum of the binary at 
the "effective" ISCO and, as such, it cannot be dependent on 
the initial separation. 

A very important consequence of our assumptions is that 
afi n for a black-hole binary is already fully determined by the 
set of coefficients S4, S5, to, t 2 , t$ computed to derive expres- 
sion (f2]). The latter, in fact, is simply the final spin for a special 
set of values for the cosine angles; since the fitting coefficients 
are constant, they must hold also for generic binaries. 

In view of this, all that is needed is to measure \£\ in terms 
of the fitting coefficients computed in refs. ifTTifOl . This can 
be done by matching expression ( TTOb with (0 [with the condi- 
tion ©J for parallel and aligned spins (a = j3 = 7 = 0), for 
parallel and antialigned spins (a = 0, j3 = 7 = tt), and for an- 
tiparallel spins which are aligned or antialigned (a = j3 = it, 
7 = or a = 7 = tt, j3 = 0). This matching is not unique, 
but the degeneracy can be broken by exploiting assumption 
(iv) and by requiring that \£\ depends linearly on cos a, cos (3 
and cos 7. We therefore obtain 
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We now consider some limits of expressions ([8]) and dTTb . 
First of all, when q — * 0, © and (ITTb yield the correct EMRL, 
i.e., |a nn | = |ai|. Secondly, for equal-mass binaries having 
spins that are equal and antiparallel, ([8]) and ( ITTb reduce to 
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Po ~ 0.687. 



(12) 



This result allows us now to qualify more precisely a com- 
ment made before: because for equal-mass black holes which 
are either nonspinning or have equal and opposite spins, the 
vector \£\ does not depend on the initial spins, expression ( TT2b 
states that \£\ M f n /4 = \£\M 2 /4 = \£\M X M 2 is, for such sys- 
tems, the orbital angular momentum at the effective ISCO. We 
can take this a step further and conjecture that \£\M\M 2 = \£\ 
is the series expansion of the dimensionless orbital angular 
momentum at the ISCO also for unequal-mass binaries which 
are either nonspinning or with equal and opposite spins. The 
zeroth-order term of this series (namely, the term 2 V3MiM 2 ) 
is exactly the one predicted from the EMRL. We note that 
although numerical simulations do not reveal the presence 
of an ISCO, the concept of an effective ISCO can neverthe- 
less be useful for the construction of gravitational-wave tem- 
plates jum. 

Finally, we consider the case of equal, parallel and 
aligned/antialigned spins (|a 2 | = |oi|, a = 0, /3 = 7 = 0, tt), 
for which expressions dTOb and dTTb become 

afin = |ai| cos/3 [1 + 1/(34 |ai| cos/3 + t + s 5 i/)] + 

v{2V3 + t 2 v + t 3 is 2 ) , (13) 
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Table I: Initial parameters of the new misaligned AEI binaries. 

where cos/3 = ±1 for aligned/antialigned spins. As expected, 
expression dT3l l coincides with (|2]i when |ai| cos/3 = a and 
with ([TJi [through the coefficients ©J when g = 1 and 
2|<2i|cos/3 = oi + 02- Similarly, ( fTOb and ( ITTb reduce 
to (O for equal, antiparallel and aligned/antialigned spins 
(|a 2 | = |ai|, a = 0, /3 = 0, 7 = 7r, or /3 = 7r, 7 = 0). 

The only way to assess the validity of expressions ([8]l 
and dTTb is to compare their predictions with the numerical- 
relativity data. This is done in Figs. Q] and [2] which col- 
lect all of the published data, together with the three addi- 
tional binaries computed with the CCATIE code |9] and re- 
ported in Table Q] In these plots, the "binary order num- 
ber" is just a dummy index labelling the different config- 
urations. The left panel of Fig. Q] in particular, shows 
the rescaled residual, i.e., (|a n n|fit — | afin | num.) X 100, for 
aligned binaries. The plot shows the numerical-relativity data 
with circles referrin g to equ al-mass, equal-spin binaries from 
refs. OIEIEISISIII Jriangles to equal -mass, unequal- 
spin binaries from refs. Illiflal. an d squares to unequal-mass, 
equal-spin binaries from refs. Il3l \TA [lit EH]. Although the 
data is from simulations with different truncation errors, the 
residuals are all very small and with a scatter of ~ 1%. 

A more stringent test is shown in the right panel of Fig.Q] 
which refers to misaligned binaries. In the top part, hexagons 



Fig ure 1 : Left panel: Rescaled residual for aligned binaries. The circles refer to equal-mass, equal-spin binaries presented in refs. JTTL[l2l . [T3l 
1 1 7L fTsL fTgll . triangles to equal-mass, unequal-spin binaries presented in ref. UHQjl], and squares to unequal-mass, equal-spin binaries presented 
in refs. Il3l.[l^. ll8l[l9ll . Here and in the right panel the "binary order number" is just a dummy index labelling the different configurations. 
Right panel: The top part reports with asterisks the final spin computed for misaligned binaries. Hexagons refer to data from ref. Jw (labelled 
"RIT"), squares to the data Table|l](labelled "AEI"), circles to data from ref. lf20h (labelled "FAU"), and triangles to data from ref. 12111 (labelled 
"PSU-UTA"). Note that these latter data points refer to the aligned component a| n since this is the only component available from ref. l2lll . 
The bottom part of this panel shows instead the rescaled residuals for these misaligned binaries. 

indicate the numerical values for |a nn | from ref. JUl, squares 
the ones in Table [fl circles those from ref. I20I1 and triangles 
those from ref. 12111 : note that these latter data points refer to 
the aligned component a| n since this is the only component 
available from ref. I2H1 . The agreement is again very good, 
with errors of a couple of percent (see bottom part of the same 
panel), even if the binaries are generic and for some the initial 
and final spins differ by almost 180° H]. 

Finally, Fig. [2] reports the angle between the final spin vec- 
tor and the initial orbital angular momentum 6* nn using the 
same data (and convention for the symbols) as in the right 
panel of Fig. Q] Measuring the final angle accurately is not 
trivial, particularly due to the fact that the numerical evolu- 
tions start at a finite separation which does not account for 
earlier evolution of the orbital angular momentum vector. The 
values reported in JH] (and the relative error-bars) are shown 
with hexagons, while the squares refer to the binaries in Ta- 
ble U and have been computed using a new approach for the 
calculation of the Ricci scalar on the apparent horizon l24ll . 
Shown with asterisks and circles are instead the values pre- 
dicted for the numerical data (as taken from refs. IH l20ll2lll 
and from Table Q} by our analytic fit (asterisks) and by the 
point-particle approach suggested in ref. Ht] (circles). 

Clearly, when a comparison with numerical data is pos- 
sible, the estimates of our fit are in reasonable agreement 
with the data and yield residuals in the final angle (i.e., 
(^fin)fit — (fen)num.) which are generally smaller than those 
obtained with the point-particle approach of ref. JH]. How- 
ever, for two of the three binaries from ref. iH the estimates 
are slightly outside the error-bars. Note that the reported an- 
gles are relative to the orbital plane at a small initial binary- 
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Figure 2: Using the same data (and convention for the symbols) as 
in the right panel of Fig. Q] we here report the angle between the fi- 
nal spin vector and the initial orbital angular momentum 8^. Shown 
instead with asterisks and circles are the values predicted for the nu- 
merical data (as taken from refs. I8l [20ll2lll and from Tabled by our 
analytic fit (asterisks) and by the point-particle approach suggested 
in ref. U (circles). 

separation, and thus are likely to be underestimates as they 
do not take into account the evolution from asymptotic dis- 
tances; work is in progress to clarify this. When the compari- 
son with the numerical data is not possible because 8^ is not 
reported (as for the data in ref. Il2l|p . our approach and the one 



in ref. J5I yield very similar estimates. 

In summary: we have considered the spin vector of the 
black hole produced by a black-hole binary merger as the sum 
of the two initial spins and of a third vector, parallel to the ini- 
tial orbital angular momentum, whose norm depends only on 
the initial spin vectors and mass ratio, and measures the orbital 
angular momentum not radiated. Without additional fits than 
those already available to model aligned/antialigned binaries, 
we have measured the unknown vector and derived a formula 
that accounts therefore for all of the 7 parameters describing 
a black-hole binary inspiralling in quasi-circular orbits. The 
equations ([8} and ( fTTT i. encapsulate the near-zone physics to 
provide a convenient, but also robust and accurate over a wide 
range of parameters, determination of the merger product of 
rather generic black-hole binaries. 

Testing the formula against all of the available numerical 
data has revealed differences between the predicted and the 
simulated values of a few percent at most. Our approach is 
intrinsically approximate and it has been validated on a small 
set of configurations, but it can be improved, for instance: by 
reducing the \ 2 of the fitting coefficients as new simulations 
are carried out; by using fitting functions that are of higher- 
order than those in expressions (fl]i and (|2}; by estimating J rad 
through PN expressions or by measuring it via numerical sim- 
ulations. 
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